Ionic correlations at the nanoscale: inversion of selectivity in a bio-nanochannel. 
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Here we show, combining a simulation and theoretical study, that electrostatic correlations typical of multi- 
valent ions can reverse the selectivity of a biological nanochannel. Our results provide a physical mechanism 
for a new, experimentally observed phenomenon, namely the inversion of the selectivity of a bacterial porin (the 
E. Coli OmpF) in presence of divalent and trivalent cations. Also, the differences and similarities between the 
driving force for this phenomenon and other similar nano and micro-escale electrokinetic effects (e.g. inversion 
of streaming current in silica nanochannels) are explored. 
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In recent years, the development of statistical-mechanical 
theories of ionic correlations near strongly charged interfaces 
0JH1 allowed a unified physical understanding of new phe- 
nomena arising in electrokinetic processes as diverse as ionic 
transport in nanochannels (3H3, DNA condensation |8[[9) or 
particle electrophoresis (TOj [TT1- Theories typically focus on 
the high electrostatic coupling regime, appearing in the case 
of multivalent ions and highly charged interfaces [T}|31. In 
this regime, ion-ion correlations are largely independent of the 
chemical nature of the interface so the physics of the problem 
can be understood using simple models for the interface and 
primitive models for the electrolyte 1 1 1. A different theoretical 
proposal shows the possibility of strong electrostatic correla- 
tions between multivalent ions and charged chemical groups 
located at the interface lfT2HT4l . These correlations are an 
interfacial analogue of the well-known Bjerrum correlations 
appearing between multivalent ions in bulk electrolyte! 12]. 
They are relevant for interfaces containing (or covered by) 
well- separated charged interfacial groups which interact in- 
dividually with ions in solution (see Figure [T]). It is important 
to note that these ion-interfacial group correlations do not re- 
quire a highly charged interface and they depend strongly on 
the nature of the structure and charge distribution of the chem- 
ical groups present at the interface. 

In this letter, we will show how correlations of the type de- 
scribed in Fig. [lj> are behind a new experimental phenomena, 
namely the inversion of selectivity observed in a biological 
nanochannel! 131 [161 in presence of multivalent cations. This 
biological nanochannel is the so-called OmpF ionic channel 
found in the outer membrane of E. coli. This is a relatively 
wide ionic channel in the sense that it allows the simulta- 
neous permeation of both cations and anions (in hydrated 
form) at high rates. This channel is made of three identi- 
cal monomerical nanopores (see Figure [2]), each one leav- 
ing an hourglass-shaped aqueous pore with a diameter be- 
tween 1-4 nm surrounded by many charged titratable residues 
with an overall negative charge 1 17] at pH=7. In presence 
of monovalent electrolyte, it has a slight cationic selectivity 
(i.e. the flux of cations is larger than the flux of anions across 
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FIG. 1: Ionic correlations near charged interfaces, (a) High electro- 
static coupling 1 1 - 3 1 . Counterions of charge +Ze accumulate near 
a strongly charged interface within a thin layer of typical thickness 
\i — kBTe/27r(—a)Ze of the order of the ionic size. At this layer, 
the cations have strong positional correlations, leaving correlation 
holes of radius ttR 2 ~ Ze/(—a) (R ^> \i so the coupling param- 
eter is large F = R/2fi ^> 1 ) (b) Correlations between ions and 
a surface fuctionalyzed with charged groups fl3lfl4l . Here, corre- 
lations result from the electrostatic interactions between multivalent 
positive counterions (green) and negatively charged oxygen atoms 
(red) located in interfacial groups. In this case, the direct interaction 
between charged interfacial atoms and counterions is essential. 



the pore) which is well understood from basic electrostatic 
concepts JT7l [T8). Unexpected experimental results |[T5l [T6l 
show that the cationic selectivity of the channel found for 
1:1 electrolytes is reversed for sufficiently high concentra- 
tions of 2:1 or 3:1 electrolytes. The different divalent cations 
(Mg 2+ , Ca 2+ , Ba 2+ , Ni 2+ ) give almost identical results, 
whereas trivalent cations (La 3+ ) induce a substantially larger 
effect. Another interesting piece of evidence is provided by 
X-ray| 19 ] data obtained from crystalls of OmpF protein chan- 
nels in 1 M of MgCl2. In this structure, hydrated Mg 2+ 
ions were found in contact with certain negatively charged 
residues, suggesting a strong cation-interfacial group inter- 
action. All this evidence (strong dependence on valency but 
small dependence on the particular cation, strong role of par- 
ticular chemical groups) suggests the existence of important 
ion-interfacial group correlations of electrostatic origin in this 
system. 

In order to identify the basic physical mechanism under- 
lying the observed selectivity inversion, we have performed 
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FIG. 2: Snapshots of the simulated system from different views, (a) top view of the trimeric OmpF protein inserted in a lipid bilayer. The 
structure of the protein has been constructed from the X-ray crystallographic data [25 1 and using the protonation states described by Varma (26l 
for the titrable residues (the resulting overall charge is —lie per monomer) (b) Front view of a snapshot showing water and ions (MgCb in 
this case). The spacing between the two water slabs corresponds to the lipid membrane and the protein (not shown here for clarity), (c) Zoom 
showing one of the channel pores and ions located inside the pore and in nearby solution (for clarity of the representation, all water molecules, 
many lipids and the other two monomers are not shown). The picture has been produced using VMD|28|. 



TABLE I: Ionic fluxes (number of ions observed to cross the OmpF 
channel during production runs of duration £ run ) and occupancy num- 
bers (defined as number of ions of each species inside the protein 
channel averaged over the simulation runs, statistical errors in the 
mean are estimated from 2a). The results correspond to simulations 
in three different conditions: 1 M KC1, 1 M MgCb and a mixture of 
1 M MgCl 2 and 1 M KC1. 





£run 


Ionic Flux 


Channel Occupancy 




(ns) 


CI" K + Mg 2+ 


CT K + Mg 2+ 


KC1 


24.9 


38 47 


21.0 ±0.2 35.5 ±0.2 


MgCl 2 


36.8 


32 - 1 


47.1 ±0.9 - 33.9 ±0.3 


Mixture 


31.6 


35 14 7 


24.4 ±0.2 7.5 ±0.2 7.2 ± 0.2 



all-atom molecular dynamics (MD) simulations with NAMD2 
l20l of the OmpF channel in different electrolytes (see Fig[2]). 
In order to study the transport properties of the channel, we 
have applied an external electric field of 14.22 mV/nm per- 
pendicular to the lipid bilayer, which creates an electrostatic 
potential drop of « 200 mV across the membrane+channel 
system (see Supporting Online material). This value of the 
potential drop was selected in order to be high enough to en- 
sure ion permeation across the channel during the simulation 
and also because it is experimentally accessible. It has to 
be stressed that full-atomistic MD simulations of transport in 
protein channels are extremely challenging, becoming possi- 
ble only due to recent improvements in algorithms and com- 
puter power (23l [241 • In this case, we have conducted the 
first simulations of ionic transport in a protein channel in di- 
valent electrolyte, which extraordinarily increases the need for 
longer simulation times and larger simulation systems, as dis- 
cussed in detail in Ref [21 ]. All simulation details, algorithms, 
processing of the data and simulation movies can be found in 
the online supplementary information l22l . 

First of all, let us note that our results for the transport of 



KC1 across the channel (see Table [I]) are in agreement with 
previous work lfTTl l24l [26l[27ll . The flux of K+ is larger than 
that of Cl~ and the average number of K + inside the channel 
is larger than the average number of Cl~ (Table)!]), as expected 
for a negatively charged channel. 

The situation is completely different in the case of MgCl2 
(see Table [I]). The channel has inverted its selectivity: a chan- 
nel conducting a slightly cationic current (in the case of KC1) 
has turned into an almost purely anion conducting channel in 
the case of MgCl2 . Let us remark that, as expected for a neg- 
atively charged channel, there is an excess of cations in its 
interior as compared with bulk electrolyte (we observe a ratio 
of 0.72 between the average numbers of Mg 2+ and Cl~ ions 
compared with 0.5 in bulk electrolyte). In spite of the excess 
of cationic charge inside the channel, its contribution to the 
current is extremely low and the vast majority of the observed 
current is due to Cl~ . We also recall that in simulations (2U 
of 1 M bulk solutions of MgCl2, 43% of the current is due to 
Mg 2+ , so the effect observed here has to be attributed to the 
protein channel. The lack of flow of the many Mg 2+ cations 
populating the channel is due to the strong, attractive inter- 
actions with channel walls. The strong correlations between 
the positions of Mg 2+ cations and negatively charged acidic 
residues of the channel can be clearly appreciated in Figure 
[3] Mg 2+ dwell near the strong concentration peaks observed 
in Fig [3^ and they are hardly found in the rest of the channel. 
In fact, residence times of Mg 2+ ions near negatively charged 
groups could be very large (of the order of 10 ns), as shown 
in the examples in Figure |3j3. In table [TT] we show the aver- 
age number of Mg 2+ in close proximity of negatively charged 
acidic groups located at the protein walls. Note that some of 
these acidic groups have an average number of Mg 2+ larger 
than 0.5, indicating that the cationic charge overcompensates 
the -e charge of the acid. Therefore, we observe local charge 
inversion [13|. 
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FIG. 3: Interaction between Mg 2+ cations and the protein channel in simulations with 1 M MgCk. (a) Concentration of Mg 2+ inside one 
of the monomers of the protein channel. We show three different isosurfaces averaged over a 10 ns trajectory fragment of the production run. 
Strong inhomogeneities (with concentrations 1 order of magnitude larger than bulk concentration) are found near negative groups located at 
the protein surface, schematically indicated as black balls (picture produced with VMD |28|) (b) Examples of long residence times of Mg 2+ 
ions near negatively charged residues evidenced by the time evolution of the distance between three selected Mg 2+ ions and certain protein 
residues indicated in the onsets. Top panel: ion close to the residue Glull7 of monomer P2 at the start of the production run which remains 
there roughly 10 ns. Medium panel: ion initially at bulk solution enters inside the P2 monomer of the channel and remains in contact with the 
Aspll3 residue during 9 ns. Bottom panel: ion initially at bulk solution enters inside the P3 monomer of the channel and remains in contact 
with both the Aspl 13 and Glu 117 residues (no detachment was observed). 



The strong correlations found between Mg 2+ ions and 
acidic residues located at the channel's wall can be consid- 
ered the first clear evidence for the transversal correlations 
predicted by recent statistical-mechanical theories of ionic 
correlations 1 1 2-H41 . In these theories, the electrostatic inter- 
action between an ion of charge qc and an interfacial atom 
of charge qi generates an excess concentration of the ion near 
the interfacial charge characterized by a pairing (or binding) 
constant Kj given by: 



Kt ~2tt 



J dn 



r ^ e QiQch 



dr. 



(1) 



In Eq.{T|, Ib = e 2 /47re/c#T is the so-called Bjerrum length 
(0.714 nm for water at 25C), do is the distance of clos- 
est approach between ions and interfacial charged atoms and 
dc = qcQilB/2 is the typical correlation length between the 
ions and the interfacial charges. In our simulations, the in- 
terfacial atoms are oxygen atoms from acidic residues (which 
have typical partial charges qc = — 0.7e). For Mg 2+ , we ob- 
tain d c « 0.5 nm. In our simulations, Mg 2+ is found to retain 
their hydration water and a distance of closest approach about 
do ~ 0-4 nm is observed, so we obtain Kj « 0.7 M _1 . In 
the case of K + , we obtain d c « 0.25 nm which is smaller 
than the sum of the crystallographic radius of oxygen and K + 
(0.28 nm), hence correlations are negligible and Kj vanishes. 
It is also interesting to estimate Kj for trivalent cations. In 
the case of La 3+ , we have d c « 0.75 nm. Since this ion has 



a size close to Mg 2+ and also tends to remain hydrated, we 
can also take d c « 0.4 nm obtaining a very high affinity of 
Kj « 6.4M -1 . Experimental results! 16] show that the in- 
version of selectivity for OmpF is larger in presence of LaCl3 
than in the MgCl2 case, as expected from our calculations. 

Further insight on the physical basis of this inversion phe- 
nomenon can be obtained by analyzing the effect of adding a 
high amount of screening monovalent salt to a situation with 
inverted selectivity. The two different kinds of ionic corre- 
lations illustrated in Figure [T] behave in opposite ways under 
the addition of high amounts of background monovalent elec- 
trolyte. Experiments in silica nanochannels found selectivity 
inversion (3 in presence of MgCi2 but the original selectiv- 
ity of the channel is recovered|5 ] (no inversion) in mixtures 
of multivalent electrolyte and KC1. This finding is in agree- 
ment with theoretical predictions! 5 1 based on ion-ion corre- 
lations between multivalent ions near strongly charged inter- 
faces (see Figure [T^). On the other hand, theories predict that 
the short range electrostatic correlations between multivalent 
ions and interfacial groups (see Figure [TJ)) are essentially not 
affected by the addition of 1 : 1 electrolyte. In this case, one 
expects that the original selectivity is not recovered after the 
addition of 1 : 1 salt; rather one may expect to observe a weak- 
ened inverted selectivity. Actually, this is what we observe 
in our simulations of the ion channel bathed by a mixture of 
1M KC1 and 1 M MgCl 2 (see Table [I]). The strong transver- 
sal correlations between Mg 2+ cations and acidic groups are 



TABLE II: Average number of Mg 2+ found nearly in contact to certain acidic residues (see Figure 3), evaluated by averaging the number 
of Mg 2+ cations at distances between do w 0.4 nm and d c w 0.5 nm of oxygen atoms of the acidic residues (see the main text) during all 
production runs. We show the results for simulations with 1 M of MgCb and a mixture of 1 M of MgCb and 1 M KC1 (statistical errors in the 
mean are estimated from 2a.) 





Glu2 


Glul81 


Asp221 


Asp74 


Aspl21 


Aspll3 


Glull7 


MgCl 2 
Mixture 


0.51 ±0.05 
0.44 ± 0.06 


0.71 ±0.05 
0.48 ± 0.06 


1.04 ±0.07 
0.65 ± 0.06 


0.28 ±0.03 
0.31 ±0.06 


0.66 ±0.05 
0.74 ± 0.08 


0.40 ±0.05 
0.50 ±0.06 


0.33 ±0.03 
0.44 ±0.05 



also observed, with a similar intensity to that observed in the 
previous case with only MgCl2 (see table [TT] and the isoden- 
sity figure available as Online Supporting Material). The ob- 
served anionic current is larger than the cationic current: the 
initial cationic selectivity of the channel is not recovered. The 
cationic current obtained with 1 M of MgCl2 and 1 M of KC1 
is significantly smaller than that obtained in presence of only 
1 M of KC1. Also, both the occupancy number and flux of 
Cl~ observed in this case are close to the values observed in 
the case with only KC1, in spite of having 3 times more an- 
ions in the electrolyte solution (see table [I]). Hence, we can 
say that selectivity effects are weakened by the addition of 1 : 1 
salt (thus demonstrating their electrostatic origin) but the in- 
version of selectivity has not been suppressed by adding 1:1 
salt. This different response to added 1 : 1 electrolyte between 
our system and silica nanochannels emphasizes the different 
nature of the ionic correlations found in each case. 

Overall, our calculations provide a strong evidence show- 
ing the role played by electrostatic correlations in the selec- 
tivity of ionic transport in a wide biological nanochannel. In 
the view of recent nanotechnological applications of biologi- 
cal nanochannels and their engineered versions (modified by 
directed mutagenesis) ll29l it is tempting to suggest that our 
results may have also implications in the design of responsive 
channels for biotechnological applications. From the point 
of view of biophysics, it is also interesting to note that recent 
electrophoretic results in phospholipid vesicles [ 30 ] also show 
a complex behavior in presence of multivalent ions which is 
robust against addition of 1 : 1 electrolyte. All these evidences 
suggest that biologically relevant molecules (lipids or pro- 
teins) have mechanisms to preserve a strong interaction with 
multivalent ions even in presence of physiologically amounts 
of screening salt. 
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